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This  paper  describes  a  novel  method  to  measure  directly  effective  diffusivity  in  electrodes  as  a  function  of 
temperature  and  relative  humidity  (RH)  at  conditions  that  are  relevant  for  proton  exchange  membrane 
fuel  cells  (PEMFCs).  The  efficacy  of  this  method  to  measure  effective  oxygen  diffusivity  (D^)  is  demon¬ 
strated  with  measurements  of  a  series  of  electrodes  of  varying  the  ionomer-to-carbon  weight  ratio  (I/C 
ratio).  The  measured  decreases  sharply  with  increasing  I/C  ratio  from  0.5  to  1.5  at  the  same  RH, 

and  reduces  gradually  with  increasing  RH  from  0%  to  100%  at  the  same  I/C  ratio.  The  measured  is 
considerably  smaller  than  the  calculated  one  using  the  Bruggeman  correction,  indicating  the  Bruggeman 
correction  drastically  underestimates  the  tortuosity  with  increasing  I/C  ratio  in  PEMFC  electrodes. 

©  2009  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

The  heart  of  a  PEMFC  is  a  three-layer  assembly  comprised  of 
a  proton-conductive  membrane  sandwiched  between  anode  and 
cathode  electrodes.  Electrodes  typically  contain  Pt  nanoparticles  as 
catalyst,  carbon  black  as  support  for  Pt  nanoparticles,  and  ionomer 
as  proton  (H+)  conductor  and  binder.  Electrodes  must  be  suffi¬ 
ciently  porous  to  transport  gases  and  water  vapor  to  carry  out 
electrochemical  reactions  inside  them.  To  improve  the  perfor¬ 
mance  of  PEMFCs,  it  is  essential  to  identify  and  understand  different 
voltage  losses.  The  cell  voltage  {EceU)  of  PEMFCs  can  be  represented 
by  the  following  equation  [1  ]: 

Ecell  —  Erev  ~  *7hOR  —  l??ORRl  —  1  '  (Ee-  +  +  EH+  An  +  EH+  Ca) 

—  tJtx(gas)’  (1) 

where  Erev  is  the  thermodynamic  potential,  ijuor  is  the  anode  over¬ 
potential  due  to  hydrogen  oxidation  reaction  (HOR),  rj0RR  is  the 
cathode  overpotential  due  to  oxygen  reduction  reaction  (ORR), 
i  is  the  applied  current  density,  Re-  is  the  electronic  resistance, 
RH+  mem *  ^h+  An *  and  ^h+  Ca  are  the  proton  resistances  in  the  mem¬ 
brane,  anode,  and  cathode,  respectively,  r]tx^gas^  is  the  gas-diffusion 
overpotential  caused  by  the  hydrogen  and  oxygen  concentration 
gradients  in  the  electrodes  due  to  gas  transport  resistances.  Except 
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for  r]tX(gCLS),  these  voltage-loss  terms  can  be  adequately  determined 
based  on  various  electrochemical  measurements  [2-5].  r]tx^gas^  is 
often  split  into  two  components:  a  dry  component  that  depends 
on  the  microstructure  of  the  porous  gas-diffusion  media  (GDL)  and 
electrode  layers  that  the  gas  must  transport  through,  and  a  wet 
component  that  is  a  function  of  the  operating  condition  combined 
with  the  ability  of  the  cell’s  components  to  reject  liquid  water. 
Clearly,  the  effective  diffusivity  of  the  electrode  layer  and  its  depen¬ 
dence  on  RH  are  important  to  fully  characterize  the  gas  transport 
resistance  of  PEMFCs. 

To  our  knowledge,  there  is  no  reported  method  capable  of 
directly  measuring  gas  effective  diffusivity  of  PEMFC  electrodes 
under  the  conditions  relevant  to  the  PEMFC  operating  conditions. 
The  Bruggeman  correction  which  is  commonly  used  to  estimate  gas 
effective  diffusivity  is  reported  to  likely  lead  to  overestimation  for 
electrodes  with  low  porosities  [6,7].  While  one  can  use  traditional 
analysis  techniques  such  as  mercury  intrusion  porosimetry  (MIP) 
and  nitrogen  adsorption  with  Barret-Joiner-Helenda  (BJH)  method 
to  characterize  the  porosity  of  electrodes,  these  techniques  are 
restricted  to  dry  analysis  conditions.  Given  that  the  ionomer  in  the 
electrode  is  highly  hygroscopic  and  swells  upon  humidification,  it 
is  important  to  be  able  to  characterize  the  electrode  layer’s  porosity 
under  conditions  more  relevant  to  fuel  cell  operation.  Additionally, 
the  thickness  and  loading  of  the  electrode  layer  (typically  lOpum 
thick  and  0.4  mgPt  cm-2,  respectively)  combine  to  make  it  difficult 
to  use  representative  samples  in  traditional  analysis  techniques. 
The  novel  method  described  herein  addresses  these  shortcomings 
by  measuring  the  effective  diffusivity  in  an  electrode  layer  under 
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Nomenclature 

List  of  parameters 

Qiir 

the  molar  oxygen  concentration  in  the  air  channel 
(mol  cm-3) 

cx 

the  molar  oxygen  concentration  in  the  nitrogen 
channel  (mol cm-3) 

Q=o 

the  molar  oxygen  concentration  in  the  inlet  of  the 
nitrogen  channel  (mol cm-3) 

Cx-L 

the  molar  oxygen  concentration  at  the  exhaust  of 
the  nitrogen  channel  (mol cm-3) 

Cy 

the  molar  oxygen  concentration  inside  the  electrode 
in  the y  direction  as  indicated  in  Fig.  3  (mol cm-3) 

Dt 

the  effective  oxygen  diffusivity  (cm2  s-1 ) 

Do2 

the  oxygen  free  space  diffusivity  (cm2  s-1 ) 

L 

the  length  of  the  nitrogen  channel  (19  cm) 

m 

the  diffusive  length  between  the  air  and  nitrogen 
channels  (0.2  cm) 

Jy 

the  molar  diffusion  flux  of  oxygen  from  both  air 
channels  to  the  nitrogen  channel  (mol  cm-2  s-1 ) 

a 

the  gas  volumetric  flow  rate  in  the  nitrogen  channel 
(cm3  s-1) 

Ax 

the  length  of  the  illustrated  segment  in  Fig.  3  (cm) 

8 

the  thickness  of  the  electrode  (cm) 

£ 

the  void  fraction  of  the  electrode 

X 

the  tortuosity  of  the  electrode 

Lpt 

the  Pt  loading  in  the  electrode  (mgPt  cm"2ectrode) 

^carbon 

the  carbon  loading  in  the  electrode 

(mgCarbonCm^ctrode) 

Monomer 

the  ionomer  loading  in  the  electrode 

(mSionomer  cmelectrode^ 

PPt 

Pt  density  (21.5 gem-3) 

Pcarbon 

carbon  density  (2.2  g  cm-3 ) 

Pionomer 

ionomer  density  (2.0  g  cm-3 ) 

conditions  similar  to  those  of  an  operating  PEMFC  and  at  a  sample 
scale  of  50  cm2.  The  efficacy  of  the  developed  method  is  demon¬ 
strated  in  measurements  of  oxygen  effective  diffusivity  (D^)  in 
electrodes  with  a  range  of  ionomer-to-carbon  weight  ratios  (I/C 
ratio). 

2.  Measurement  of  inside  electrodes 

°2 

2.1.  Method  description 

The  measurement  is  based  on  oxygen  “in-plane”  diffusion  in 
electrodes  on  a  conventional  PEMFC  platform  with  50  cm2  flow 
fields.  Fig.  1  shows  the  schematic  cross-section  of  the  cell  fixture  for 
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Fig.  1.  Schematic  cross-section  of  the  cell  fixture. 
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Fig.  2.  Plan-view  of  the  flow  field. 

the  measurements.  The  “in-plane”  diffusion  means  that  the  dif¬ 
fusion  direction  is  parallel  to  the  electrode  surface  (the  x  direction 
in  Fig.  1 ).  In  operating  PEMFCs,  oxygen  is  more  likely  to  experience 
“though-plane”  diffusion  in  electrodes,  where  the  diffusion  direc¬ 
tion  is  perpendicular  to  the  electrode  plane  (they  direction  in  Fig.  1 ). 
However,  in  porous  PEMFC  electrodes  should  be  the  same  no 
matter  oxygen  is  in  either  “in-plane”  or  “through-plane”  diffusion. 
Fig.  1  shows  the  porous  electrode  film  coated  on  ethylene  tetraflu- 
oroethylene  (ETFE)  is  tested  by  mounting  it  between  a  flow  field 
with  three  parallel  gas  channels  and  a  blank  flow  field.  Air  flows 
through  the  two  outer  channels  and  nitrogen  through  the  center 
channel.  This  arrangement  results  in  oxygen  diffusion  from  the  air 
channel  into  the  nitrogen  channel  via  the  porous  electrode.  Because 
the  nitrogen  channel  is  flanked  on  both  sides  by  the  air  channels, 
the  rate  of  oxygen  diffusion  is  effectively  doubled,  thereby  reduc¬ 
ing  the  required  length  of  flow  channel  to  achieve  a  measurable 
oxygen  concentration  at  the  exhaust  of  the  nitrogen  channel. 

A  gasket  made  of  silicone  rubber  with  the  thickness  of  250  p,m 
(Durometer  35A,  McMaster-Carr)  is  used  to  seal  the  air  and  nitrogen 
channels  under  a  load  of  5.3  kN.  The  gasket  is  cut  by  laser  to  match 
the  channels.  The  width  of  the  gasket  that  seals  the  lands  between 
the  air  and  nitrogen  channels  are  2  mm  as  shown  in  Fig.  1 ,  allowing 
some  tolerance  for  placement  of  the  gasket  to  ensure  that  it  is  com¬ 
pletely  supported  by  the  3  mm  wide  flow  field  lands.  Placement  of 
the  gasket  is  important  because  the  width  of  the  gasket  between 
the  air  and  nitrogen  channels  defines  the  oxygen  diffusion  length 
used  in  the  calculation  of  De^. 

Fig.  2  shows  the  plan-view  of  the  flow  field  depicting  the  fea¬ 
tures  of  the  air  and  nitrogen  channels.  Air  and  nitrogen  co-flow 
through  the  flow  field  with  a  channel  length  of  19  cm.  The  oxygen 
concentration  at  the  exhaust  of  the  nitrogen  channel  is  monitored 
by  a  self-heated  oxygen  sensor  (model:  LZA03-E1 ,  NGK  Spark  Plugs, 
Michigan).  The  test  is  run  with  the  two  gas  streams  held  at  ambi¬ 
ent  pressure  ( 1 01 .3  kPa).  The  pressure  drops  in  the  air  and  nitrogen 
channels  are  calculated  to  be  21  and  17  Pa,  respectively.  The  max¬ 
imum  pressure  gradient  between  the  air  and  nitrogen  channel  is 
only  4  Pa  when  both  air  and  nitrogen  are  flowing.  Even  if  a  pressure 
gradient  of  1 0  kPa  exists  between  the  air  and  nitrogen  channels,  the 
resulting  oxygen  bulk  convection  is  calculated  to  account  for  less 


Z.  Yu,  R.N.  Carter  /  Journal  of  Power  Sources  195  (2010)  1079-1084 


1081 


o 

II 

X 

x  =  L  (19  cm) 

\  z 

Air-  -  Cjw* - ► 

Jv 

*  —  2 

r__  x+  Ax 

XT  Cx  1 

N2  — ->j 

c 

[ _ 

1 

i  t  !/«  =  2mm 

_ \L _ k _ 

Air - - ► 


Fig.  3.  Illustration  of  the  oxygen  mass  balance  over  the  nitrogen  channel. 


than  2.1 3%  of  the  oxygen  diffusion  flux  in  the  measurements  for 
electrodes  with  I/C  ratios  of  0.5  and  1.0.  Thus,  the  driving  force  for 
forced  convective  flow  between  the  adjacent  air  and  nitrogen  chan¬ 
nels  is  negligible,  and  the  transport  of  oxygen  from  the  air  channels 
to  the  nitrogen  channel  is  driven  only  by  the  oxygen  concentration 
gradient.  Because  of  the  minimal  pressure  gradient  between  the 
external  air  channel  and  the  ambient,  it  is  assumed  that  there  is  no 
much  oxygen  in  the  air  channel  diffusing  through  the  electrode  to 
the  ambient.  Thus,  oxygen  can  only  diffuse  from  the  air  channels  to 
the  nitrogen  channel  though  the  porous  electrode. 

To  measure  the  oxygen  diffusion  length  accurately  under  the 
load,  a  piece  of  Fuji  Prescale  film  for  super  low  pressure  (FujiFilm, 
Japan)  is  inserted  between  the  electrode  and  the  gasket  and  then 
compressed  by  the  two  flow  fields.  Under  the  compression,  the  Fuji 
Prescale  film  develops  colorful  strips  due  to  the  pressure  exerted 
by  the  gasket  between  the  air  and  nitrogen  channels.  The  width  of 
colorful  strips  on  the  pressure-paper  reflects  that  the  expansion  of 
the  2  mm  wide  gasket  under  a  load  of  5.3  kN  is  negligible.  Thus,  the 
oxygen  diffusion  length  is  2  mm  in  the  calculation  of  D^. 


Because  the  air  and  nitrogen  channels  are  kept  at  the  same  tem¬ 
perature,  pressure,  and  RH,  and  oxygen  bulk  convection  between 
the  air  channels  and  the  nitrogen  channel  is  negligible  as  discussed 
above,  oxygen  transport  from  the  air  channels  to  the  nitrogen  is 
assumed  to  be  an  oxygen-nitrogen  binary  diffusion.  Thus  Jy  can  be 
expressed  according  to  Fields  law  of  diffusion  as  follows: 


I  _  _neff  ~~~y 
Jy  ~  °2  dy  ’ 


dCv 


(4) 


where  is  the  effective  oxygen  diffusivity,  Cy  is  the  molar  oxy¬ 
gen  concentration  inside  the  electrode  in  the  y  direction  in  Fig.  3. 
In  the  present  study,  the  test  conditions  and  fixture  design  are 
chosen  such  that  the  measurement  of  PEMFC  electrodes  would 
result  in  oxygen  concentrations  of  less  than  400  ppm  in  the  nitro¬ 
gen  channel,  which  is  500  times  lower  than  21%  oxygen  in  the  air 
channels.  Thus,  the  oxygen  concentration  in  the  air  channel  (Cair) 
is  considered  to  be  constant  at  21%  although  a  trace  amount  of 
oxygen  diffuses  from  the  air  channels  to  the  nitrogen  channel.  Con¬ 
sequently  ,Jy  can  be  assumed  to  be  constant  down  the  length  of  the 
channels  because  the  oxygen  concentration  gradient  (i.e.,  21%  vs. 
at  most  400  ppm)  between  the  air  channels  and  the  nitrogen  chan¬ 
nel  is  basically  unchanged.  Thus,  Jy  in  Eq.  (4)  can  be  simplified  as 
follows: 


Jy  =  -D 


eff  Coir 
°2  m  ’ 


(5) 


where  Cair  is  the  molar  oxygen  concentration  in  the  air  channel 
(defined  by  the  operating  temperature  and  RFI)  and  m  is  the  diffu¬ 
sive  length  scale  defined  by  the  sealing  gasket  width  (2  mm). 

Thus,  with  a  constant Jy  and  a  boundary  condition  as  Cx= o  =  0,  the 
oxygen  concentration  at  the  exhaust  of  the  nitrogen  channel  (Cx=/J 
can  be  derived  from  Eq.  (3)  as  follows: 


(6) 


Then,  combination  of  Eqs.  (5)  and  (6)  yields  the  following  expres¬ 
sion  for  calculating  : 


2.2.  Calculation  ofD\ ^ 

As  described  above,  oxygen  is  designed  to  “in-plane”  diffuse 
through  a  2  mm  wide  porous  electrode  due  to  the  oxygen  concen¬ 
tration  gradient  between  the  air  channels  and  the  nitrogen  channel. 
Since  the  electrode  thickness  ( 1 0  p,m)  is  200  times  smaller  than  the 
oxygen  diffusion  length  (2  mm),  the  oxygen  concentration  inside 
the  electrode  is  assumed  to  be  the  same  in  the  y  direction  in  Fig.  1. 
In  addition,  the  oxygen  concentration  inside  the  electrode  which  is 
under  the  air  or  nitrogen  channel  is  assumed  to  be  the  same  in  the 
x  direction  in  Fig.  1 .  Therefore,  the  oxygen  transport  in  this  study  is 
simplified  to  an  issue  of  oxygen  diffusion  due  to  the  concentration 
gradient  through  a  porous  electrode  which  is  2  mm  wide,  10p,m 
thick,  and  19  cm  long. 

Fig.  3  illustrates  an  oxygen  mass  balance  over  a  segment  of  the 
nitrogen  channel  for  developing  the  equation  in  the  calculation  of 
Dq^.  It  is  expressed  as  follows: 

0.  x  (Q+Ax  -  Cx)  +  2Jy  x  Ax  x  <5  =  0  (2) 

where  Q.  is  the  gas  volumetric  flow  rate,  Cx  is  the  molar  oxygen 
concentration  in  the  nitrogen  channel ,Jy  is  the  molar  diffusion  flux 
of  oxygen  from  both  air  channels  to  the  nitrogen  channel,  Ax  is  the 
length  of  the  illustrated  segment  in  Fig.  3,  and  8  is  the  thickness  of 
the  electrode.  In  the  limit  of  Ax^  0,  Eq.  (2)  becomes: 

dCx__2Jy8 

dx  a  •  1  j 


r} eff  _  Cx=lQJTI 
°2  2 8LCair  ' 


(7) 


Note  that  Cx=i ,  Cair,  and  Q  are  all  RH-dependant  variables  in  the 
calculation. 


3.  Experimental  details 

Carbon-supported  Pt  catalyst  (50  wt.%  Pt/carbon;  Tanaka,  Japan) 
and  ionomer  solutions  (5  wt.%  Nafion®  in  organic  solvents, 
900  equivalent  weight  (EW);  DuPont,  USA)  were  ballmilled  in  a 
mixture  of  water  and  ethanol  for  72  h.  The  resulting  catalyst  inks 
were  then  coated  by  Meyer  rod  coating  on  blank  ethylene  tetrafluo- 
roethylene  (ETFE)  films  at  a  catalyst  loading  of  0.4  mgPt  cm-2  (note: 
cm2  refers  to  the  geometric  area  of  electrodes).  Electrodes  with 
I/C  ratios  of  0.5, 1.0,  and  1.5  were  fabricated  and  then  die-cut  into 
50  cm2  samples. 

The  electrodes  are  mounted  in  the  cell  fixture  shown  in  Fig.  1 
and  then  compressed  at  a  load  of  5.3  kN,  giving  1.41  MPa  pressure 
on  the  porous  electrode  film.  The  cell  fixture  is  heated  to  80  °C.  Air 
and  nitrogen  flow  rates  into  the  saturators  are  300  and  150  seem, 
respectively.  Air  and  nitrogen  are  vented  to  ambient  after  the  flow 
field.  The  connection  between  saturators  and  the  flow  field  as  well 
as  between  the  flow  field  and  the  oxygen  sensor  are  heated  to 
prevent  water  condensation  in  those  regions.  This  is  critical  for 
maintaining  the  desired  RH  in  the  air  and  nitrogen  channels  as  well 
as  preventing  artifacts  in  the  apparent  oxygen  diffusion  due  to  the 
presence  of  liquid  water  in  the  cell  fixture.  The  temperature  of  the 
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Fig.  4.  Results  of  the  control  test  showing  the  oxygen  concentration  in  the  exhaust 
of  the  nitrogen  channel  as  a  function  of  compressive  load  applied  to  the  cell  fixture. 
The  test  sample  is  a  piece  of  uncoated  ETFE.  Test  conditions  are  80  °C,  300/1 50  seem 
for  air/nitrogen,  101  kPa,  and  0%  RH. 


Fig.  5.  The  last  two  cycles  of  oxygen  concentrations  at  the  exhaust  of  the  nitro¬ 
gen  channel  (C*=i)  with  RH  sweep  from  20%  to  100%  over  an  electrode  with  I/C  of 
0.5.  Ten  cycles  of  RH  sweep  are  conducted  over  each  electode  to  get  stable  oxygen 
concentrations. 


saturators  is  adjusted  to  control  the  RH  in  air  and  nitrogen  chan¬ 
nels  to  values  of  20%,  40%,  60%,  80%,  90%,  and  100%.  Additionally, 
the  gas  streams  can  be  bypassed  around  their  respective  satura¬ 
tors  to  obtain  0%  RH.  Ten  cycles  of  RH  sweep  from  20%  to  100%  are 
conducted  over  each  electrode  to  get  stable  measurements.  Fol¬ 
lowing  the  ten  cycles  of  RH  sweep,  the  measurement  at  0%  RH  is 
carried  out  for  5  h  to  get  stable  value.  The  oxygen  concentration  in 
the  exhaust  of  the  nitrogen  channel  is  recorded  every  10  s.  At  each 
condition,  the  oxygen  concentration  is  averaged  over  stable  data  in 
a  period  of  30  min.  Three  repeats  with  new  samples  are  conducted 
to  determine  the  repeatability  of  this  method.  Scanning  electron 
microscopy  (SEM,  FEI  Quanta  400)  is  used  to  measure  the  thickness 
of  electrodes  before  and  after  the  tests.  The  fresh  electrodes  on  ETFE 
are  embedded  in  a  standard  epoxy  resin  (Epoes,  Struers)  for  the 
cross-section  SEM  image,  while  the  tested  electrodes  transferred 
onto  silicone  gasket  are  self-supported  for  the  same  analysis. 

4.  Results  and  discussion 

A  gasket  is  used  to  seal  the  air  and  nitrogen  channels.  A  con¬ 
trol  test  was  performed  to  determine  what  compressive  load  is 
necessary  to  achieve  adequate  sealing  without  inducing  excessive 
stress  on  the  porous  electrode  structure  and  potentially  causing 
deformation.  In  the  control  measurement,  an  ETFE  film  without  the 
electrode  coating  is  mounted  in  the  cell  fixture  to  test  the  extent  of 
sealing  with  increasing  compression  load  at  80  °C,  101  kPa  and  0% 
RH.  Fig.  4  shows  the  results  of  this  test,  where  the  oxygen  concen¬ 
tration  in  the  exhaust  of  the  nitrogen  channel  is  plotted  as  a  function 
of  the  compression  load.  The  oxygen  concentration  decreases  from 
about  18  to  10  ppm  as  the  compression  load  increases  up  to  5.3  kN, 
where  upon  it  remains  fairly  steady  up  to  8.0  kN.  Considering  the 
oxygen  diffusivity  in  silicone  rubber  is  typically  in  the  order  of 
10-7  to  10-6cm2  s-1  [8],  permeation  of  oxygen  through  the  sili¬ 
cone  gasket  should  result  in  an  oxygen  concentration  of  less  than 
0.2  ppm  in  the  nitrogen  channel.  This  analysis  indicates  that  the 
background  oxygen  concentration  of  1 0  ppm  is  mostly  due  to  small 
leaks  at  the  gasket  seal.  The  oxygen  concentration  in  the  nitrogen 
exhaust  in  the  tests  for  typical  electrodes  (i.e.,  I/C  ratio  from  0.5  to 
1.0)  ranges  generally  from  120  to  350  ppm,  which  is  many  times 
larger  than  the  1 0  ppm  baseline.  Thus,  a  compression  load  of  5.3  kN 
on  flow  fields  is  chosen  in  this  method  to  achieve  adequate  seal¬ 
ing  for  the  air  and  nitrogen  channels.  The  oxygen  concentration  of 
10  ppm  is  assumed  as  a  baseline  and  subtracted  from  the  oxygen 
concentrations  measured  in  the  electrode  tests. 

Fig.  5  presents  the  last  two  cycles  of  oxygen  concentrations  at  the 
exhaust  of  the  nitrogen  channel  {Cx=l  in  ppm)  at  varying  RH  values 


over  an  electrode  with  I/C  ratio  of  0.5.  This  result  demonstrates  that 
Cx=l  is  a  steady-state  value  at  each  RH  condition  and  drastically 
sensitive  to  the  RH  change  over  the  electrode.  Note  that  CX=L  needs 
to  be  converted  into  the  unit  of  molar  concentration  for  calculating 
in  Eq.  (7),  thereby  incorporating  the  dilution  effect  of  increasing 
RH  into  the  calculation. 

Fig.  6  displays  SEM  images  of  the  cross-section  of  electrodes  with 
an  I/C  of  1.0  before  and  after  the  measurement.  The  fractures 
observed  in  the  tested  electrode  are  caused  by  sample  preparation 
for  SEM  analysis.  The  SEM  images  indicate  no  change  in  the  elec¬ 
trode  thickness,  thereby  showing  the  porous  electrode  does  not 
collapse  under  the  5.3 1<N  compression  load.  More  SEM  analysis  (not 
presented)  shows  that  electrodes  with  I/C  ratios  of  0.5  and  1.5  also 
have  the  same  thickness  of  10  pm  after  the  measurements. 

Thus,  a  thickness  of  10  pm  is  used  in  Eq.  (7)  to  calculate  for  all 
the  three  types  of  electrodes. 

Fig.  7  displays  the  measured  for  electrodes  with  I/C  ratios 
of  0.5,  1.0,  and  1.5  as  a  function  of  RH.  The  results  show  that  the 
measured  of  electrodes  decreases  sharply  with  increasing  I/C 
ratio  at  the  same  RH  but  reduces  gradually  with  increasing  RH 
at  the  same  I/C  ratio.  For  example,  at  0%  RH,  the  measured 
decreases  70%  and  95%  as  I/C  ratios  increases  from  0.5  to  1 .0  and  1 .5, 
respectively.  In  contrast,  the  measured  of  the  electrodes  with 
I/C  ratios  of  0.5  and  1.0  correspondingly  decrease  20%  and  25%  as 
RH  increases  from  0%  to  100%.  Given  that  the  electrode  thickness  is 
observed  in  SEM  analysis  to  be  independent  of  the  I/C  ratio,  increas¬ 
ing  the  I/C  ratio  from  0.5  to  1 .5  triples  the  ionomer  volume  fraction 
and  correspondingly  reduces  the  electrode  void  fraction.  Due  to  the 
ionomer’ s  nature  of  swelling  upon  hydration,  one  expects  to  see  a 
decrease  in  the  diffusivity  as  RH  is  increased.  The  measured  of 
the  electrode  with  I/C  ratio  of  1 .5  is  significantly  smaller  than  those 
of  electrodes  with  I/C  ratios  of  0.5  and  1.0,  indicating  the  severe 
ionomer  filling  and  blocking  of  pores  inside  the  electrode  with  I/C 
ratio  of  1.5  for  oxygen  diffusion.  Furthermore,  the  measured 
of  the  electrode  with  I/C  ratio  of  1.5  appears  independent  on  RH 
because  (i)  its  absolute  value  is  too  small  to  show  the  dependence 
and  (ii)  the  developed  method  may  lose  its  sensitivity  when  mea¬ 
suring  the  value  of  CX=L  in  the  range  of  15-32  ppm  with  a  10  ppm 
baseline  for  the  oxygen  concentration. 

Fig.  7  also  presents  the  remarkable  reproducibility  among  the 
three  replicate  tests  for  each  type  of  electrodes,  further  demon¬ 
strating  the  efficacy  of  the  developed  method  for  measuring 
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Fig.  6.  SEM  images  of  the  cross-section  of  (a)  a  fresh  electrode  with  I/C  of  1 .0  coated 
on  ETFE  and  (b)  a  tested  electrode  with  I/C  of  1.0  supported  on  the  silicone  rubber 
gasket.  Note  that  the  tested  electrode  film  always  transfers  to  the  silicone  gasket 
after  being  compressed  in  the  cell  fixture. 


in  PEMFC  electrodes.  Stumper  reported  that  in  combination 
of  the  gas-diffusion  layer  and  the  catalyst  layer  was  measured  to 
be  between  0.0061  and  0.0072  cm2  s-1  [9].  Although  he  did  not 
report  the  catalyst  loading  or  I/C  ratio  in  the  catalyst  layer,  his  mea¬ 
surements  are  reasonably  close  to  the  results  repeated  here  for  an 
electrode  with  I/C  ratio  of  1.0. 


Fig.  7.  Measured  De^  (average  of  three  replicates)  as  a  function  of  RH  for  elec¬ 
trodes  with  I/C  ratios  of  0.5,  1.0,  and  1.5.  Test  conditions  are  80  °C,  300/1 50  seem 
for  air/nitrogen,  and  101  kPa.  The  error  bars  indicate  95%  confidence  limits. 


Table  1 

Summary  of  measured  ,  calculated  porosities,  and  calculated  D^a. 


Electrodes 

Measured  D^fb 
°2 

Calculated  porosity 

Calculated  Djf 

°2 

I/C  0.5 

0.02410 

68.1% 

0.155 

I/C  1.0 

0.00734 

58.1% 

0.122 

I/C  1.5 

0.00119 

48.1% 

0.092 

a  Conditions:  80  °C,  101  kPa,  and  0%  RH. 
b  Unit:  cm2  s_1. 


According  to  its  definition,  is  expressed  as  follows: 

<  =  fo02,  (8) 

where  s  is  the  void  fraction,  r  is  the  tortuosity,  and  Dq2  is  the  oxygen 
free  space  diffusivity.  Typically,  the  Bruggeman  expression  is  used 
to  estimate  the  tortuosity  as  [7]: 

r  =  g-0-5.  (9) 

Combination  of  Eqs.  (8)  and  (9)  give  the  Bruggeman  correction  for 

De£  in  PEMFC  electrodes: 

°2 

Di(  =s'-5Do2.  (10) 

In  the  02-N2-H20  ternary  system  in  PEMFC  electrodes,  Dq2  is  usu¬ 
ally  estimated  by  the  Wilke  equation  [10-12],  with  the  assumption 
of  zero  flux  of  nitrogen  and  water  vapor  in  the  system: 

U  _ _ 1  ~  X°2 _  (^^) 

°2  (Xn2/Do2-N2)  +  (xH2o/^02-H2o)’ 

where  x02 ,  xN2 ,  and  xH2o  are  the  mole  fractions  of  oxygen,  nitrogen, 
and  water  vapor  in  electrodes,  respectively.  The  binary  diffusivity  of 
oxygen  and  nitrogen,  Dq2  _n2,  is  calculated  using  Chapman-Enskog 
equation  [10]  and  the  binary  diffusivity  of  oxygen  and  water  vapor, 
0o2-h2o.  is  calculated  using  Slattery-Bird  equation  [13].  At  80 °C, 
101  kPa,  and  0%  RH,  Dq2  is  equivalent  to  Dq2_ n2  and  calculated  to 
be  0.276  cm2  s-1. 

The  porosity  of  electrodes  at  0%  RH  can  be  determined,  since  the 
electrode  thickness  as  well  as  the  mass  loadings  and  densities  of  Pt, 
carbon,  and  ionomer  are  known.  The  expression  for  the  calculated 
porosity  of  electrodes  is 

g  _  ^  _  (^Pt/ PPt)  +  (^carbon/ Pcarbon)  T  ( L  ionomer/  Pionomer)  (12) 

where  e  is  the  porosity  of  electrodes,  Lpt  is  the  Pt  load¬ 
ing  (0.4  mg  cmd2ectrode  in  this  study),  Lcarb0n  is  the  carbon  load¬ 
ing  (0.4  mg  cm^ctrode  for  50  wt.%  Pt/ carbon),  and  Lionomer  is  the 
ionomer  loading,  which  equals  to  Lcarbonx\IC.  ppt  (21.5 gem-3) 
[14],  Pearbon  (2.2  gem  3 )  [14],  and  Pionomer  (2  §  cm  3 )  [15]  are  den¬ 
sities  of  Pt,  carbon,  and  ionomer,  respectively.  V  is  the  geometric 
volume  of  electrodes  (0.001  cm3  cm-2ctrode  in  this  study). 

Table  1  summaries  the  measured  D^,  the  calculated  porosity, 
and  the  calculated  using  the  Bruggeman  correction  for  elec¬ 
trodes  at  80  °C,  101  kPa,  and  0%  RH.  Although  the  measured 
decreases  by  95%  as  I/C  ratio  increases  from  0.5  to  1 .5,  the  calculated 
porosities  only  reduces  by  29%.  Thus,  the  disproportionate  decrease 
in  the  measured  is  most  likely  attributable  to  an  increase  in  the 
tortuosity  in  electrodes  with  increasing  I/C  ratio.  In  addition,  the 
measured  is  considerably  smaller  than  the  calculated  for  all 

the  three  electrodes.  The  ratio  of  the  measured  DeJf  to  the  calculated 

°2 

decreases  as  I/C  ratio  increases,  indicating  the  Bruggeman  cor¬ 
rection  drastically  underestimates  the  tortuosity  with  increasing 
I/C  ratio  in  PEMFC  electrodes  [7]. 


1084 


Z.  Yu,  R.N.  Carter  /  Journal  of  Power  Sources  195  (2010)  1079-1084 


5.  Conclusion 

A  method  has  been  developed  to  measure  directly  gas  effective 
diffusivity  in  electrodes  as  a  function  of  temperature  and  RH  under 
the  conditions  relevant  to  the  PEMFC  operating  conditions.  At  80  °C 
and  101  kPa,  in  electrodes  with  I/C  ratios  of  0.5,  1.0,  and  1.5 
are  measured  on  a  50  cm2  flow  field  for  PEMFCs,  with  RH  varying 
from  0%  to  1 00%.  The  measured  of  electrodes  decreases  sharply 
with  increasing  I/C  ratio  at  the  same  RH  but  reduces  gradually  with 
increasing  RH  at  the  same  I/C  ratio.  The  calculated  porosity  of  elec¬ 
trodes  reduces  by  29%  as  I/C  ratio  increases  from  0.5  to  1.5,  but  the 
measured  decreases  by  95%,  inferring  a  dramatic  increase  in  the 

tortuosity  with  increasing  I/C  ratio.  In  addition,  the  measured  is 

considerably  smaller  than  the  calculated  using  the  Bruggeman 
correction  for  all  the  three  electrodes,  indicating  that  the  Brugge¬ 
man  correction  drastically  underestimates  the  tortuosity  in  PEMFC 
electrodes. 
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